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Abstract. We present optical i?-band light curves of the double gravitationally lensed quasar FBQ 0951+2635 
from observations obtained at the Nordic Optical Telescope between March 1999 and June 2001. A time delay of 
At — 16 ±2 days (1<t) is determined from the light curves. New constraints on the lensing geometry are provided 
by the position and ellipticity of the lensing galaxy. For a (fi m , = (0.3,0.7) cosmology, the time delay 
yields a Hubble parameter of Ho = 6Qt% (random, la) ±2 (systematic) kms -1 Mpc -1 for a singular isothermal 
ellipsoid model and Ho = 63l^ (random, la) ±1 (systematic) kms -1 Mpc -1 for a constant mass-to-light ratio 
model. In both models, the errors are mainly due to the time-delay uncertainties. Non-parametric models yield 
-ffo =64^7 (random, la) ±14 (systematic) kms -1 Mpc -1 . 
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1. Introduction 

Refsdal (:1962I) was the first to point out the possibility 
of determining Ho through gravitational lensing of a vari- 
able source. Since the light travel times for the various 
images are unequal, intrinsic variations of the source will 
be observed at different times in the images. In particular, 
the time delay between images is a measurable parame- 
ter related to the gravitational potential of the lens and 
Ho. There is also a weaker dependence on cosmological 
parameters and on the source (z s ) and lens (z^) redshifts. 
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* Based on observations made with the Nordic Optical 
Telescope, operated on the island of La Palma jointly by 
Denmark, Finland, Iceland, Norway, and Sweden, in the 
Spanish Observatorio del Roque de los Muchachos of the 
Instituto de Astrofisica de Canarias. Based on observations 
made with ESO Telescopes at the Paranal Observatories un- 
der programme ID 66.A-0062(A). Based on observations made 
with the NASA/ESA Hubble Space Telescope, obtained from 
the Data Archive at the Space Telescope Science Institute, 
which is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS 5- 
26555. These observations are associated with program #7887. 



Prompted by the successful optical measurements of 
time delays in QSO 0957+561 (Schild & Thomson 115571 
Kundic et al. IT557|) and PG 1115+080 (Schechter et al. 
1997}, we conducted photometric monitoring campaigns 
of gravitationally lensed quasars at the Nordic Optical 
Telescope (NOT) in La Palma and at the 1.5-m Danish 
Telescope in La Silla between 1997 and 2001, yielding four 
time delays. The first one was obtained for the doubly im- 
aged quasar B1600+434 which is lensed by an edge-on 
spiral (Burud et al. EUOUj) . Then Hjorth et al. ffiKESj l esti- 
mated the delay for the two (summed) components of RX 
J0911+0550. Finally, delay measurements were obtained 
for the doubly lensed quasars HE 2149—2745 (Burud et 
al. llWOikJl and SBS 1520+530 (Burud et al. 12002 bjl . These 
observations increased the total number of firm time- 
delay measurements in gravitational lenses to ten (e.g. 
Kochanek & Schechter |2TjU4"J> . 

In this paper we present the fifth result from our pro- 
gramme: the time-delay measurement in the doubly-image 
quasar FBQ 0951+2635. It was discovered by Schechter 
et al. I|1998l hereafter S98) as a double quasar with an- 
gular separation of l'.'l. The images are fairly bright with 
mean i?-band magnitudes of 16.9 and 18.0 for A and B, 
respectively. The time delay between the two images is 
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Fig. 1. A finding chart for FBQ 0951+2635. The star used 
to model the point spread function (PSF) is indicated. 
The two reference stars (labeled S2 and S3) used for the 
photometry are also marked. 

expected to be relatively short, of the order of one month. 
The monitoring programme, the light curves and the time 
delay are discussed in Sects. 00] Spectroscopic observa- 
tions with the Very Large Telescope (VLT), including a 
measurement of the source redshift and the detection of 
microlensing, are reported in Sect. 03 In Sect. we analyse 
an HST/NICMOS image to constrain a mass model of the 
galaxy and use it to estimate Ho. Finally, the main results 
are summarised in Sect. 

2. Observations and data reduction 

The system was observed on a regular basis at the NOT, 
starting in March 1999. For the first 18 months it was 
observed roughly twice per month but thereafter around 
once per week. The target is visible from the beginning 
of October to the end of June, which results in gaps of 
about ^3 months in the light curves. There were addi- 
tional gaps due to periods of bad weather and bad see- 
ing. Two different instruments were used: Andalucia Faint 
Object Spectrograph and Camera (ALFOSC) and the 
Standby CCD Camera (StanCam), equipped with detec- 
tors yielding pixel scales of 0'.'188 and 0'.'176, respectively. 
In Fig.Qthe central region of the field of view of ALFOSC 
is displayed. One data point typically consisted of 1-3 ex- 
posures of five minutes each. The seeing varied from 0'.'7 to 
2'.' 8 during the roughly two years of observations, with 0'.'9 
being the most frequent value. All the imaging data were 
pre-processed (bias-subtracted and flat-field corrected) us- 
ing standard IRAF routines. 

3. Photometry 

The photometric data consist of one stacked frame per 
epoch. Accurate photometry of the blended quasar im- 
ages is crucial for time-delay measurements and all data 
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Fig. 2. Left: stacked i?-band image of a 12" x 12" region 
centered on FBQ 0951+2635. The seeing is ~l'/2 and the 
total exposure time is ^8.3 hours. Right: deconvolved im- 
age (FWHM = 0"28) obtained from the simultaneous de- 
convolution of 49 frames. North is up and east to the left. 

were analysed with the versatile deconvolution algorithm 
developed by Magain, Courbin & Sohy l|1998l hereafter 
MCS). This algorithm has already been used to analyse 
the data of several blended quasar images (e.g. Burud et 
al. ESTinni l2TY)2ll I2002bl Hjorth et al. EISSJ- The main ad- 
vantage of using the MCS deconvolution in the case of a 
photometric monitoring of a target is its ability to decon- 
volve all the frames from different epochs simultaneously. 
This constrains the positions of the two quasar images, 
which do not vary with time. The intensities, however, 
are allowed to vary from one image to the other, hence 
producing the light curves. 

The two quasar components are clearly distinguished 
in our deconvolved image of FBQ 0951+2635 (see Fig.0), 
but the lensing galaxy remains too faint to be detected. 
Although this is clearly a disadvantage when measuring 
the redshift of the galaxy (Sect. 03), it is fortuitous for the 
time-delay estimate since a contamination by an extended 
object could complicate the analysis. 

The light curves of FBQ 0951+2635 consist of 58 data 
points in the i?-band, as presented in Fig. 03 The magni- 
tudes are calculated relative to the magnitude of the point 
spread function (PSF) star (R — 16.6) shown in Fig. ^ 
Two other stars were deconvolved as well to check for sys- 
tematic errors. The error bars include both photon noise 
and additional systematic errors, e.g. PSF errors. The lat- 
ter are estimated by using the reference stars, as explained 
in Burud et al. ijMTnjl . 

4. Time-delay estimate 

It is clear from Fig. that the quasar did not show large 
intrinsic variations. But a direct comparison of the first 
five data points of the two components indicates that the 
time delay is very short since the gradual rise in the B 
component closely follows the rise in the A component. In 
addition, during the last ^100 days the light curves display 
a nearly simultaneous decrease in flux. 

During the observations, a possible tiny short-lived fea- 
ture was detected in both components as shown in the top 
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Fig. 3. The A and B component i?-band light curves for 
FBQ 0951+2635 and two reference stars in the field. The 
plotted values are R-0.9 for B, R- 1.8 for S2 and R-2.0 
for S3. The magnitudes are calculated relative to the cal- 
ibrated PSF star in the field. The error bars represent 
photon noise and PSF errors estimated from the decon- 
volution of the two reference stars. Since S2 and S3 are 
significantly fainter than components A and B, their pho- 
ton noise is larger. The dashed horizontal lines represent 
the mean magnitudes of stars S2 and S3. 



panel of Fig.0J Sliding the light curves across one another 
to find the best visual agreement, and making use of this 
feature, lead us to a rough estimate for the time delay 
of ^15 days. The light curves show that A is the leading 
component. Once light curves of quasar components are 
obtained, statistical methods are needed to determine the 
time delay more accurately. Using the \ 2 minimization 
method described in Burud et al. (|2D01 ) , a more objective 
value of At = 17 ± 4 days (la) is found from the light 
curves. The magnitude difference between the A and B 
images is 1.01 ± 0.01 mag, corresponding to a flux ratio 
between the A and B components of Fa/Fb — 2.54 + 0.02. 
The errors quoted here are obtained from Monte Carlo 
simulations of 1000 sets of light curves, assuming that the 
photometric errors are uncorrelated and follow a Gaussian 
distribution. 

The method used to reach the aforementioned esti- 
mate adds a linear term to one of the components, thus 
modelling slow microlensing effects (in Sect. we find in- 
dependent evidence for microlensing in FBQ 0951+2635 
by analysing its spectra). This linear correction does not 
remove all the external variations as can be seen in the 
middle panel of Fig. 21 Faster variations on time scales of 
50-100 days are still present. The best method for dealing 
with such high-order variations is the iterative version of 
the algorithm (Burud et al. I2001|l . It was also applied to 
the light curves in an attempt to correct for these fast vari- 
ations. This method yields At = 13 ±4 days (la), slightly 
lower than the value found with the direct method but 
still in agreement within the error bars. 
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Fig. 4. Top: a zoom in on the light curves produced from 
the latest observations. The plotted values are R — 0.9 
for the B component in order to ease the comparison. A 
short-lived dip in the light curves is indicated by the two 
circles. Middle: time-delay shifted light curves of FBQ 
0951+2635. B is shifted by +1.01 magnitudes and for- 
ward in time by 17 days. A linear correction for the exter- 
nal variations is applied. Bottom: time-delay shifted light 
curves of FBQ 0951+2635. B is shifted by +1.06 mag- 
nitudes and forward in time by 16 days. Only the last 38 
data points were used with a direct method (no correction 
for external variations). 
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Table 1. Estimated time delays and i?-band magnitude 
differences for FBQ 0951+2635 calculated with three dif- 
ferent methods. The quoted uncertainties are la errors. 





At 


Am 


Method 


(days) 


(mag) 


x 2 fit 


17 + 4 


1.01+0.01 


Iterative fit 


13 + 4 




Last 38 points 


16 + 2 


1.06 + 0.01 



The uncertainties remain large because our method 
actually finds two possible time delays for the whole data 
set when applying only the linear term. In addition to the 
17 days, we also find that 25 days are consistent with the 
light curves. The result is sensitive to a smoothing parame- 
ter in the program, which smooths out dates without data 
points on the model curve. Since there are large gaps in our 
light curves, we explore various values for the parameter 
and find the latter result is less robust, with 17 days a more 
stable result. To verify this, we also ran the programme on 
the last 38 data points spanning a 260 day interval (shown 
in the top panel of Fig.0J, which do not include any large 
gaps and are more frequently sampled . In addition, this 
interval includes the potential short-lived feature. This ap- 
proach yields a time delay of At = 16 + 2 days (la), with 
no sign of a second possible minimum. The magnitude dif- 
ference between the A and B images is 1.06+0.01 mag, cor- 
responding to Fa/Fb = 2.65 ± 0.02. The resulting shifted 
light curves are displayed in the bottom panel of Fig. 0] 
In the same figure there are signs of possible microlens- 
ing variations in one of the components for the first few 
hundred days. This microlensing is most likely the result 
of compact objects in the lensing galaxy (positioned close 
to the B component; see below), located in the line-of- 
sight towards the quasar images. This result is further 
strengthened by the independent spectroscopic microlens- 
ing observations in Sect. [5] 

The time-delay estimates and magnitude offsets ob- 
tained from our three different approaches are consistent 
with each other and are presented in Tabled Given that 
the large gaps in our light curves affect our results, we 
adopt At = 16 + 2 days as the best estimate for the time 
delay. 

5. VLT spectroscopy 

The redshift of the lensing galaxy has proved elusive. S98 
found absorption lines attributed to Mgn at z — 0.73 
and z = 0.89 in both spectra of the quasar components. 
From the position of the lens on the fundamental plane, 
Kochanek et al. ( 2000) suggested 2d = 0.21 with a range 
of 0.18-0.23, while Kochanek et al. ifTMfl hereafter K98) 
quote a redshift of 0.24, estimated from the properties of 
the lensing galaxy (C. Kochanek 2002, private commu- 
nication). Multiply lensed quasars have also become an 
important tool to detect microlensing by compact objects 



within the lensing galaxy. This offers the opportunity to 
study the spatial structure of the lensed quasar with very 
high angular resolution. 

On 2001 March 24-25 we obtained deep optical spec- 
tra of FBQ 0951+2635 with the VLT using the Multi 
Object Spectroscopy (MOS) capability of FORS1 with l'/4 
wide slitlets. Two grisms were used, G300V+GG435 cov- 
ering approximately 4300-8600 A and G300I+GG435 with 
a wavelength coverage of roug hly 6600-11000 A. 

We acquired 2 x 1200 + 620 s spectra with G300V and 
3 x 1200 s spectra with G300I. The slitlets were oriented 
with a position angle of —55? 6 to cover both quasar com- 
ponents and the lensing galaxy. One of the slitlets was cen- 
tred on a star in the field of view whose PSF was later used 
as an input for deconvolution of the quasar spectra. The 
individual spectra were combined, yielding a seeing of 0'.'4 
(0"6) and a spectral resolution of 10 A (16 A) FWHM for 
G300I (G300V). The projected pixel size was O'.'l x 2.59 A 
(O'.'l x 2.69 A) for G300I (G300V). 

Redshift estimations of lensing galaxies are often lim- 
ited by observational difficulties arising from the fact that 
the light from these galaxies can be heavily blended with 
the light from one of the quasar images, the latter some- 
times being several magnitudes brighter. The spectral de- 
convolution algorithm developed by Courbin et al. (2000) 
is able to deconvolve spectra in the spatial direction. It 
can therefore be used to separate the quasar and galaxy 
spectra. This technique was applied to FBQ 0951+2635 
in order to deblend the spectra of the two quasar images 
and to search for the lensing galaxy, located only 0'.'2 away 
from component B. 

In the spectra of both quasar components we identify 
three emission lines at >5<r, namely Cm] A1909, Mgn 
A2798 and [O n] A3727. These lines place the quasar at 
z s = 1.246 + 0.001, a value fully consistent with that found 
by White et al. (2000). In addition, the absorption lines 
first reported by S98, interpreted as the Mgn doublet at 
z = 0.729 and z — 0.892, are detected in the A component. 
In the spectrum of the B component we only detect the 
z = 0.892 system. 

Strong equivalent width (EW) difference exists be- 
tween the Mg II line in the spectra of the two quasar im- 
ages, with EW a factor of 1.46 smaller in component B 
than in component A with high significance. The most 
probable explanation is that the continuum emission of 
the B component is microlensed by compact objects in 
the lensing galaxy. This result is consistent with standard 
models of quasar structure, with the scale of the broad 
emission line region significantly greater than that of the 
continuum source (e.g. Wambsganss et al. 119901 Krolik 
EHII. The Mgn flux ratio between A and B is 3.7 ± 0.3. 
This is considerably higher than the continuum flux ratio 
(~2.5) and that from broadband photometry. 

In order to obtain the lensing galaxy spectrum, free 
of any contamination by the two bright nearby quasar 
images, the 2-D spectrum of FBQ 0951+2635 was spa- 
tially deconvolved. However, the signal-to-noise ratio was 
very low in the lensing galaxy spectrum. Due to this and 
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Fig. 5. Left: an HST/NICMOS image of FBQ 0951+2635. Right: a deconvolved HST/NICMOS image shown on a 
logarithmic intensity scale. In both panels the lensing galaxy is denoted by the letter 'G'. 



a lack of significant absorption and emission lines in the 
spectrum, we were unable to determine a spectroscopic 
redshift. In what follows we adopt z& = 0.24 from K98, 
but note that decreasing z& from 0.26 to 0.18 decreases 
the value of Ho by 37%. 

6. Lens modelling 

The time delay and redshifts measured for FBQ 
0951+2635 can be used to infer an estimate of Hq, based 
on modelling of the total gravitational potential responsi- 
ble for the lensing effect. This includes the main lensing 
galaxy and any tidal perturbations from objects near it 
or along the line of sight. FBQ 0951+2635 is not an easy 
case to model, as the lensing galaxy is not readily visible. 
However, as seen in the left panel of Fig. [SI there is a small 
'lightbump' northwest of, but close to the B component. 

We use two distinct approaches to convert the time 
delay into Hq: parametric models that involve an analyti- 
cal form for the lensing galaxy, and non-parametric mod- 
els that involve more degrees of freedom, i.e. a pixelated 
surface density (Saha & Williams 119971 Williams & Saha 
2000 , Saha & Williams EMI) . 

6.1. HST/NICMOS results 

In March 1998, four 640 s images of the FBQ 0951+2635 
lens system were obtained by K98, taken with the 
Near Infrared Camera and Multi-Object Spectrometer's 
(NICMOS) NIC2 camera onboard the HST. The filter 
used was F160W. The images were drizzled to give a pixel 
scale of 0'.'0375. 

We deconvolved the drizzled image, using the A com- 
ponent as an empirical PSF. The right panel of Fig. 
shows the result from the deconvolution: the lensing 
galaxy is positioned close to the B component. This lens 
configuration is in full agreement with the fact that A is 



Table 2. Image and lens positions. 



Object 


ft. 


0y 


Lens 


O'.'OOO ± 0'.'002 


O'.'OOO ± 0'.'002 


A 


Of.' 750 + 0'.' 001 


C'459 ± O'.'OOl 


B 


-0'.' 142 + 0'.' 001 


-0'.'169± O'.'OOl 



the leading component (see Sect.^J, as predicted by grav- 
itational lens theory. We fit ellipses to the isophotes in 
the core of the lensing galaxy resulting in an ellipticity of 
e = 0.25 ± 0.04 and a position angle (PA) of 22° ± 4°. 
Assuming that the mass distribution follows the light dis- 
tribution we can restrict ourselves to these ranges of el- 
lipticities and PAs. For the de Vaucouleurs model (see 
below) we also estimated the effective radius, resulting in 
R e = C'09 ± 0'.'02. The exact positions of the quasar im- 
ages are listed in Table [21 in a lens-centered coordinate 
system. The quasar images have an angular separation 
of 1'.'091 ± O'.'OOl, while the B component is positioned 
0'.'221 ± 0'.'002 from the center of the lensing galaxy. We 
note that our position of the galaxy obtained from the 
HST/NICMOS image is quite different from that pub- 
lished by S98, where a lens candidate is detected mid- 
way between the two components. But the S98 data were 
ground-based and they acknowledge that the candidate 
might be a spurious detection. 

6.2. Analytical models 

Using the lensmodel package (Keeton 2001), we fit the lens 
system with a dark matter model and a constant M/L 
model. The former is a singular isothermal ellipsoid (SIE) 
model with a flat rotation curve, while the latter is the pro- 
totypical constant M/L model: the de Vaucouleurs model. 
This is a similar approach as Kochanek & Schechter ( 2004 ) 
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Fig. 6. Ax 2 contours for FBQ 0951+2635 in the plane of 
the ellipticity (e) and shear (7). The contours are drawn 
at Ax 2 = 2.30, 6.17 and 11.8, corresponding to the la, 
2a and 3<r confidence intervals for two degrees of freedom. 
There is a weak degeneracy between e and 7. The area 
between the horizontal dashed lines is the range of the 
measured ellipticity. 



carried out for four lens systems with known time de- 
lays measured by Burud et al. l|2UUUl I2002al I2002bjl and 
Schechter et al. J1997). The SIE model sets a lower bound 
on H , while the constant M/L model sets an upper bound 
on Hq . This results from the fact that no galaxy in a cold 
dark matter dominated Universe should have mass distri- 
butions which are more centrally concentrated than the 
luminosity distribution (Mo et al. 199BJ). In addition, a 
typical lensing galaxy should not have a rising rotation 
curve at approximately twice the effective radius, making 
the SIE the least concentrated reasonable model. 

During the fit, the astrometry of the quasar images 
relative to the lens is fixed, as well as the redshift of the 
source and lens. The models are constrained by using the 
range of ellipticities and PAs obtained in Sect. 16. II In addi- 
tion, we use Fa/Fb as a constraint. Four values have been 
derived independently. In Sect. ^ we found from broad- 
band photometry that Fa/Fb = 2.5-2.7. This is similar to 
the ratio calculated from the spectral continua (Sect.EJ). 
On the other hand, the emission line flux ratio is measured 
to be Fa/Fb = 3.7 + 0.3 (Sect.GJ. Finally, S98 report the 
radio emission line ratio to be 4.7 ± 0.6. While flux ratios 
from broadband photometry are generally known to be 
unreliable due to the ever present possibility of microlens- 
ing, the emission line fluxes should reflect the "true" flux 
ratio (e.g. Wisotzki et al. 2003). Hence, in what follows, 
we adopt Fa/Fb = 3.7 ± 0.3 as a model constraint, but 
note that allowing Fa/Fb to vary between 2.5 and 4.7 
does change the value of Ho by only approximately 4%. 

The free parameters in our model include the mass of 
the lensing galaxy, the source coordinates and flux, and 
the amplitude of the external shear (7) along with the 
direction to it. We include shear in our models since at 
least two mass concentrations lie along the line of sight, 
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Fig. 7. Top: the slope of the projected mass density pro- 
file plotted against Hq for 121 different non-parametric 
models. Bottom: probability distribution of the derived Hq 
values for the same models as above. The median value of 
the distribution is indicated by the dotted line. 

namely the ones responsible for the Mg 11 absorption sys- 
tems (S98; Sect. |5J). It is important to understand any 
degeneracies between the two sources of angular struc- 
ture in the model, the ellipticity and shear. The external 
shear breaks the circular symmetry of the lens and there- 
fore it has the same effect as introducing ellipticity in the 
lens. Figure (BJ shows the contours of the A\ 2 in the space 
of e and 7 after optimizing all other parameters. In the 
SIE model this results in a moderate shear of 7 = 0.07, 
while for the constant M/L model there is a large shear of 
7 = 0.29. By evaluating Hq inside the la contours in Fig. El 
we get an estimate for the systematic uncertainties in both 
models. If we take the currently most popular cosmology, 
i.e. fi m = 0.3, Qa — 0.7, and At = 16 days, we find 
i?o = 60^y (random, la) ±2 (systematic) km s" 1 Mpc -1 
for the SIE model, while Ho — 63^y (random, la) ±1 (sys- 
tematic) kms -1 Mpc^ 1 for the constant M/L model. Due 
to the relatively low redshift of the lensing galaxy, the 
influence of a change in cosmology is small compared 
to the measurement errors. Adopting other cosmologies, 
(fi ro , fi A ) = (1.0, 0.0) or (O m , fi A ) = (0.3, 0.0), will change 
the inferred value of Hq by only about 4% for this lens 
system. 

6.3. Pixelated models 

While the mass distribution of the pixelated models of 
Saha & Williams l|2004[l is not as well related to the physi- 
cal parameters of the lens, the models allow an exploration 
of a wide range of lens shapes. Thus, this approach better 
constrains the systematic errors as it is unrestricted by 
the confines of the parametric models. 

We fit the lens system with the PixeLens programme. 1 
During the fit, the ellipticity of the mass distribution was 

http: / /ankh-morpork. maths. qmw.ac.uk/~saha/ astron/lens/pix/ 
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kept free. As in the analytical models, we added external 
shear. The lensing galaxy appears unperturbed in Fig. 
leading us to enforce the inversion symmetry constraint 
(i.e. the lens looks the same if rotated by 180°). In order to 
avoid models less concentrated than the SIE, we imposed 
a constraint stating that the index of the projected radial 
mass density profile (a) must be steeper than —1. As there 
is not an option in PixeLens to invoke this constraint, we 
ran a total of 500 models and, in the upper panel of Fig.[7| 
plot those 121 that fulfill a < —1. Most lens systems, 
when modelled independently, show a correlation between 
a and Hq, in the sense that steeper density profiles re- 
sult in higher values of Ho- FBQ 0951+2635 shows this 
correlation rather weakly. The lower panel of Fig. [7\ dis- 
plays the derived Hq probability distribution. The mean 
of the distribution is calculated to be Hq = 64ly (random, 
la) ±14 (systematic) km s _1 Mpc -1 . 

7. Discussion 

The time delay At = 16 ± 2 days (la) has been mea- 
sured in the lensed quasar FBQ 0951+2635 on the basis of 
i?-band images obtained with the NOT. Due to the lack 
of pronounced features in the light curves and because 
the time delay is short compared to the sampling of the 
light curves, Ar is not very well constrained, thus explain- 
ing the large error bars. If this short time delay is to be 
pinned down more accurately, intensive monitoring should 
be performed in the future. As long as the variability is 
sufficiently strong and erratic during the observations, this 
should help constrain At even further. 

Applying an SIE lens model to the case of FBQ 
0951+2635, we derive H = 60±? (random, la) ±2 (sys- 
tematic) km s _1 Mpc , where the error bars incorporate 
the uncertainty in the measured time delay as well as 
the uncertainties in the lens model due to the degener- 
acy between the shear and ellipticity. For a constant M/L 
model we estimate Ho = 63^ 7 (random, la) ±1 (sys- 
tematic) kms -1 Mpc -1 . These two analytical models were 
studied because they represent the limiting case for phys- 
ically possible mass distributions given our current un- 
derstanding of dark matter. Non-parametric models em- 
ployed on the lens system resulted in if — 64l 7 (ran- 
dom, la) ±14 (systematic) km s -1 Mpc -1 , with a rela- 
tively broad range in a. 

Unfortunately we cannot determine which mass profile 
is appropriate given the available data. However, by us- 
ing an independent Hq estimate we can try to break the 
degeneracy. The local estimate by the HST Key Project 
(Freedman et al. IM)T)> is H = 72 ± 8 km s -1 Mpc" 1 , 
the recent results by WMAP (Spergel et al. I2UU3|| indi- 
cate that Hq = 71^3 km s _1 Mpc -1 , while the value of 
H Q = 59 ±6kms -1 Mpc -1 by Saha et al. pM)l is based 
on Cepheid calibration for a sample of galaxies with Type 
la supernovae. Due to our large At uncertainty, all three 
estimates are consistent with the models we adopt. It is 
therefore important to monitor this lens system more fre- 



quently in order to beat down the random uncertainty on 
any H estimate from FBQ 0951+2635. 
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